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Abstract. The spoutaueous emergence of heterogeneous dislocatiou patterns is a 
conspicuous feature of plastic (ieforuiatioii Hud strain hardening of ciystalliiie solids. 
Despite long-standing efforts in the materials science and physics of defect couimnultles, 
there is no general consensus regarding the physical mechanism which leads to the 
formatiou of dislocation patterns. In oitler to establish the fundamental mechanism, 
v,v foimulate an extiemely simpliiled. minimal model to investigate the formation 
of patterns based ou the continuum theory of fluxes of cui^'ed dislocations. We 
demonstrate that strain hardening as embodied In a Taylor-type dislocation density 
dependence of the flow stress, In conjunction with the structure of the kinematic 
equations that govern dislocation motion under the attiou of external stresses, is 
aheady sufficient for the formation of dislocation patterns that are consistent with 
the principle of similitude. 
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1. Introduction 

Work hardening during plastic deformation of crystalliue solids is associated with 
significant changes in dislocation micrastiiictnre. The increase in dislocation density 
on the specimen scale is accompanied by the quasi spontaneous emergence of regions of 
low dislocation density and clusters of high dislocation density which to a large extent 
persist upon unloading. These metastable structures are denoted as dislocation patterns. 
Despite a significant degree of morphological variation depending on slip geometry 
and loading mode (e.g. cell [T] or labyrinth structines |2j. dislocation accumulation 
in veins |2] or walls |3|), these patterns are characterized by some fairly universal scaling 
relationships. These relationships are commonly referred to as law of similitude" or 
'similitude principle! They relate the characteristic length A of deformation-induced 
dislocation patterns to the applied stress at which they have formed, and to their 
average dislocation density p: The wavelength A is proportional to the dislocation 
spacing and inversely proportional to the applied stress, A = D/^/p = DaGb/Toxt 
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where G is the shear modulus, 6 is the modulus of the Burger's vector, a ss 0.3 is the 
uon-dimeusional coefficient relating flow stress and dislocation density in the Taylor 
relationship Tg^i = aGb^fp, and the remainiug parameter D is typically of the order of 
Z? = 10... 20. These relations have been observed to hold in a wide range of materials 
and over more than four orders of magnitude in scale |3, For a general oveiwiew, see 
e.g. the topical review by Sauzayi^Kirbin [S]. 

Already for more than half a century attempts have been made to model dislocation 
patterning, using a huge variety of different simulation methods and theoretical 
approaches (see e.g. [7] for an overview). Early models were often based on analogies 
with other physical problems including spinodal decomposition patterns [S], chemical 
patterning iu react ion-diffusion sv'stems [31 [TU|, noise-induced phase transitions [TH [TZ) . 
or even used simple internal-energy minimization arguments |ia| . Despite the large 
mmiber of published models, no consensus regarding the key meciinnisms of dislocation 
patterning has emerged. In fact, it is only too easy to criticize imicii of the published 
literature as being inconsistent with basic obser\'ations about dislocation behavior and 
dislocation patterning: Energy minimization approaches which either consider energetic 
quasi-equilibrium ("low energy dislocation structures") or patterning in the approach 
to thermal equilibriiim (spinodal decomposition) are inconsistent with the fact that 
dislocation patterning is only observed in intrinsically non-equilibrium situations where 
dislocations are driven by external stress. Analogies with reaction-diffusion models 
are spurious because it is difficiilt to see how the directed motion of dislocations in 
response to Peach-Koehler forces coiilcl be described as a diffusion process, In our view 
many of the published approaches focus too strongly on the patterning aspect, i.e., they 
represent attempts at ’modeling patterns* using templates whicli do not fully account for 
the rather peculiar dynamic properties aud interactions of dislocations. Instead, we are 
convinced that the only viable method to anive at a satisfactory theoretical description 
of dislocation patterns consists in ’modeling dislocations’, i.e., developing models that 
are capable of desenbing the stress-driven motion and interactions of curved and flexible 
dislocation lines in crystal lattice structures subjected to specifled external constraints: 
if we get the physics right, patterns are bound to emerge. 

The huge variety of modeling approaches is matched by a similarly wide range of 
opinions regarding the physical processes that drive dislocation patterning. The spinodal 
decomposition model of Holt |5| considers reduction of the elastic energy stored in the 
long-range stress fields of dislocations to be the key driver. Kulihnann-Wilsdorff and 
her school (see e.g. [T3]) more generally hold that dislocation patterns represent the 
configurations of lowest internal energy that are kinematically accessible at a given 
stage of deformation. By contrast, Mughrabi |3| points out that dislocation patterns 
evolve during deformation under the action of an applied stress and in conformity 
with constraints imposed by the deformation. In this picture, the heterogeneity in 
the dislocation distribution leads to heterogeneity in the local flow stress. The ensuing 
plastic strain heterogeneity needs to be elastically accommodated and causes substantial 
deformation-induced long-range internal stresses, accompanied by an appreciable stored 
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elastic strain energy [TJ. Finally, Madec et al. [TS| state that only short-range stresses 
(line tension), which control the foniiation of junctions, together with processes such as 
cross slip which help the entanglement of dislocations, are essential to patterning while 
long-range stress fields are irrelevant. 

To resolve these controversies it would be desirable to develop conipiitational 
models whicli allow to activate or de-activate physical micro-mechaiiLsms at will, and 
thus to decide which features of dislocation motion are essential and which ones arc 
incidental to pattern formation - an approach chosen for instance by Madec et al. 
|15| . Three-dimensional discrete dislocation dynamics (3D DDD) simulations provide 
a complete description of dislocation motion and interactions on the level of individual 
dislocation lines and might therefore offer a suitable methodological ft’nmework for sucii 
an investigation. Indeed, early stages of dislocation patterning such as the formation of 
irregular dislocation clusters, braids, or entanglements, have been frequently observed in 
simulations, see e.g. Madec et al |15| . Devincre et al. [E] and recently also Hussein et 
al. [T2|> published investigation, however, allows to clearly identify the pattern wave 
length and to directly link this through the similitude scaling relation to the applied 
stress. The reason for this is simply the required size such a simulation would need 
to have. As observed by El-Awandy et al. |18) . the deformation behavior does not 
depend on the absolute size of a [simulated or real) specimen but on the size in units of 
dislocation spacings - a fact that is immediately obvious from the scaling properties of 
dislocation sv'stems [TS|. What is a sufficient size for investigating dislocation patterns 
in a DDD simulation? Experimental dislocation patterns have typical wav'eleugths 
(spacings between cell walls or braids, veins) of the order of D = 10...20 dislocation 
spacings |5]. In order to clearly identify a pattern wavelength and to ensure tliat the 
pattern is representative of biilk behavior rather tlian governed by boundary effects, 
the system should exliibit regular features in all spatial directions, hence, its spatial 
dimension L should at least be 3-4 times larger than the pattern wavelength which makes 
for at least a linear dimension of 30-40 dislocation spacings. This is at the upper end 
of the system sizes attainable by state-of-the art 3D DDD simulations. Table 1, which 
shows values of system size and final dislocation density for some simulations reported in 
the literatme. makes this point explicit. In fact, to our knowledge, there is only a single 
example of a 3D DDD simulation in the piiblished literature where the simulated system 
size readies some tens of the final dislocation spacing [201 and even there, the achieved 
strain of ': 2% may have been insufficient to form fully developed dislocation patterns. 
Ill all other mentioned studies, tendencies towards dislocation clustering concomitant 
with the formation of dislocation-depleted regions can be observed but the scale of 
these features is comparable to the simulated volume and hence no definite conclusions 
on pattern morphology or pattern wavelength can be drawm. 

Patterning can be efficiently simulated using two-dimensional discrete dislocation 
dyuamics (2D DDD) where dislocations are envisaged as straight parallel lines. I)ut 
sucli models require phenomenological rules to represent genuinely three-dimensional 
processes associated with segment curvature, such as dislocation multiplication and 
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reference 

linear 

dimension L 

final dislocation 
density p 

effective size 

lVp 

Madec et el. |li}| 

10 ^lm 

2,0 X 10^^ 

14 

ArsenliB et al. 1501 

D pm 

7.5 X 10^^ m-2 

42 

HasKeiii et al. [T7| 

10 [iiu 

3.0 X 10^2 

17 


Table 1: Effective sizes of some DDD simulations reported in the literatiu'e 


junction formation. Attempts have been made to formulate and parameterize such 
rules using 3D simulations (’2.5D DDD'), whicli results in patterns that are consistent 
with the similitude principle |2I]. Howe^’er, it is clear that this approacli owing to its 
geometrical restrictions cannot give access to dislocation patterns in general deformation 
geometries. 

We expect that the above mentioned limitations of 3D DDD, which mostly relate 
to the huge computational cost of the method, can be overcome by resorting to 
density-based continuum dislocation dynamics (CDD) frameworks. A suitable CDD 
framework must be capable of representing the kinematics of curved dislocations in 
a geometrically correct manner and should allow to include the essential features of 
dislocation interactions. Several attempts in this direction have been reported in the 
literature. Some studies focus on subgrain formation (e.g. [22|) * a process associated 
with geometrically necessary’ dislocations only - and can therefore not capture the initial 
stages of dislocation patterning which are associated with the clustering of statistically 
stored dislocations of zero net Burgers vector. Other investigatiorrs focus on specific 
mechanisms such as the sw'eepirrg of najTow dislocation dipoles by curved dislocations 
as envisaged by Kratochvil et al. |23l which may only be relevant in specific deformation 
conditions such as fatigue. 

In the present manuscript we undertake first steps towards a generic model that 
captures the coupled dynamics of statistically stored and geometrically necessary 
dislocations while accoimting for the specific kiirematics of curved dislocation lines. To 
this errd we use the kinematic framework of Hochrainer’s CDD theory' to represent 
fluxes of cuiwed dislocations in a continuum setting [231 ESI E] We use this 
fiairrework to iirvestigate the following question: Given the kinematics, what are the 
minimal ingredients required for the emergence of dislocation patterns? In Sect. we 
formulate our kinematic CDD framework together with a minimal model of dislocation 
interactions and a specification of boundary conditions. Scaling the system reveals the 
structure of the equations and allows to analyze the critical conditions required for 
patterning in Sect. We then use numerical simulation to investigate the properties of 
the emergent patterns and their eralution with stress during strain hardening (Sect.0. 
We conclude with a discussion of the physical mechanism of dislocation patterning and 
point out perspectives tor future work. 
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2. The Continuum Model 

2.1. Continuum dislocation dynamics theory 

As kinematic framework on whicii we base our imnierical investigation we use the CDD 
theory formulated by Hochrahier and co-workers |23 EZ] (in the spirit of the present 
work one of the authors has also used the more complex, so-called higher-dimensional 
CDD' theor>' in a simple ’1,5D model' of dislocation patterning |ZS|)« The field variables 
descril)ing the dislocation microstructure on a given slip sv'stem with normal vector n 
and Burgers vector 6 = 6 s are the total density pi, the dislocation density vector q. and 
the cur\'ature density qi. The latter quantity can be envisaged as an av'eraged product 
of dislocation liue density and line cur\'ature: its integral J^q^dV = 27rN^ over a volume 
V yields the imml)er of dislocation loops Nd contained in the volume, hence, qt can also 
be envisaged as a loop density. We define the GKD density \'ector as ^ = — {l/ 6 )Vn 7 
where 7 is the plastic slip on the considered slip system, and Vn = V * (1 — n S' n) 
is the Nal)la operator in the slip plane where 1 is the rank-2 miit tensor. We assume 
that dislocation motion occurs by crystallographic glide with average (scalar) velocity 
V. Hence, Orowan’s equation |2S] gives the evohrtion of q as = ptbi\ 

We note that our definition of q differs from the definition of the dislocation density 
vector K in [!2ril[251l27] : both formulations are related by ^ = — n x k where n Is the slip 
plane nomial vector. The modified definition leads to an equivalent, but more intuitive 
formulation of the kinematic equations which are now giv'eu by 


^Ipt = 

- V„ - {vq) + vg, 

(1) 

dtQ = 

- V„(?.'p,) 

(2) 

= 

-V„-('-rQ<'>+A®-V„r). 

(3) 


These equations can be considered lowest-order terms of an alignment tensor expansion 
\27\. We close this expansion by expressing the curvature density vector and the 
rank -2 alignment tensor iu terras of the variables p^,, and g. 

The curvature density vector in (j^ plays very much the same role as the 
dislocation density vector ^ in (|^. As pointed out by Hochrainer et ai. [2U|. 
can be approximated as the product of the dislocation density vector and the mean 
cun’ature qilpi. 

Q'‘’ = -e^. (4) 

Pi 

In the vicinity of a homogeneous and isotropic dislocation airangemeut, this 
approximation becomes exact. The rank-2 aligumeut tensor has trace Tr A^^^ = pt 
and for an isotropic dislocation arrangement, A^^^ = 0,5ptl The deviatoric part Aj^?^ 
characterizes imbalances in the distribution of dislocation characters (e.g. predominance 
of edge vs. screw dislocations). For a weakly anisotropic arrangement containing 
a fraction of geometrically necessary excess dislocations, we follow the proposal of 
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Monavari et aJ. [SU]) and wiite 

where Ip is the unit v'ector in direction of the GND vector (the slip gradient) and ^ is 
a function with the properties ^(0) = 0 and ^(1) = 1. The approximation for 
proposed by Hochrainer et ai. [2S| Is tantamount to setting ^(u) = y/u which appears 
inappropriate since it leads to a formulation which is non-analytic at n = 0 , i.e., in the 
vicinity of a homogeneous and isotropic dislocation arrangement. We therefore follow' 
the suggestion of Monavari et al. |3!3 from a Maximum Entropy argument deduce 
^>{u) ft; u(l -h u)/2. 


2.2. Dislocation interactions and dislocation velocity 


ill a density-based model, dislocation interactions enter in different ways depending 
on their range. Interactions on scales above the characteristic dislocation spacing (in 
the foUow'ing referred to as /on^-range interactions) can be represented in terms of 
eigenstresses associated with a heterogeneous distribution of the plastic strain ([ 211 1 ^ ^ 
or, equivalently, with the presence of geometrically necessary dislocations (non-vanishing 
dislocation density vector g). Interactions on scales below the characteristic dislocation 
spacing, on the other hand (in the following referred to as short-range interactions) 
cannot be described in this manner, for the obvious reason that the continuous 
dislocation density fields and the associated strain field cannot meaningfully represent 
heterogeneities below the scale of the single-dislocation spacing. These short-range 
interactions are descnbed in the standard manner in terms of a friction-like yield stress, 
which we assume to depend on the local value of the dislocation density according to 
the Taylor relationship. Additional terms relating to line cun’ature or to short-range 
interactions between geometrically necessary dislocations might be considered (see the 
disciLssion of such terms by Zaiser ^ Sandfeld [T2|) but, in the spirit of formulating a 
minimal model, here we exclusively use the Taylor relationship to describe the short- 
range dislocation interactions whicli control the flow stress. The dislocation velocity is 
then computed from the glide component of the Peach-Koehler force as 


V = 


(6/B)sign(r){|T| - tj ) IF |r| > Tf, 
0 ELSE. 


( 6 ) 


where B is the drag coefficient. In this expression, r is the resolved shear stress 
in the sHp system, computed by solving the elastic boundary value problem for a 
body containing the eigenstrain field 7 , This stress captures both the effect of the 
externally applied boundary tractions and/or constraints, and the long-range dislocation 
interactions associated with a heterogeneous plastic strain or, equivalently, with the 
dislocation field g. 
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We study two vei'sions of the model, namely (1) a situation where long-range interactions 
are absent but short-range interactions are present, and (2) a situation with long-range 
interactions but no short-range ones. In both cases we envisage a layer of plastically 
active material of thickness h where deformation occurs in single slip (a slip zone or 
slip band). The layer is assumed parallel to the xy plane which is also the slip plane, 
and the slip direction is the i direction. For nimiericai simulations, periodic boundary 
conditions are assumed in the i and y directions. 

In case (1) we assume that the plastically active layer forms a thin Him loaded by 
homogeneous shear tractions acting on its unconstrained top and bottom surfaces. In 
case (2) we consider the layer to be embedded into an infinite block of elastic material 
which is loaded remotely, again in a spatially homogeneous manner. The stress field in 
the plastically active layer is then given by 


T = 7io:i, case(l), 

r = rcxL J S(i-i\y-y')-r{j:\y')dx'dy\ caae(2), (7) 


where the shear stress due to the externally applied tractions and the Fourier 

transform of the kernel 5 is given by (see e.g. |33p 




Gh \ kl 

2 


( 8 ) 


The convolution integral in (Q can from this relation be evaluated in standard manner 
by a double Fast Fourier TVansform. 


3. Mathenmaticsd analysis 
S.l. Dimensionless scaling 

To analyze the behavior of the model it is convenient to switch to a non-dimensional 
representation w'hich allows to identify the independent parameters. We define 
the scaling relations between quantities with physical \mits and their dimensionless 
counterparts (indicated by a tilde) os r = CrT (tor stresses), pt = CpPi (tor dislocation 
densities) and x = (for lengths), with the scaling factors 

a = C, = pa-. = (9) 

where po = (pi(t = 0)) is the mean initial dislocation density. The same scaling applies 
to the total dislocation density and the dislocation density vector. The scaling relation 
for the lengths implicates that lengths are measured in multiples of average dislocation 
spacings. For the velocity we find for r > Tf 
B rr 

i; = = l-TinVA- 

where Ua = oG6v^/n.xt- 


(10) 
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The scaling of velocity implies a scaling of time, since the scaling factor of velocity 
can be understood as the ratio between the characteristic length Cx a^id a characteristic 
time Cl- Flirtheimore, the fact that the dislocation curvature has the dimension of an 
inverse length suggests to scale in proportion with CpjCx. Finally, we observe that 
the plastic slip can be envisaged as the product of the Burgers vector and the area per 
unit volume swept by dislocations. These observations lead to the derivative scaling 
relations i' = Cv^s t = and 7 = C 77 with 


a = 


b 


Br 


ext 


p ^ Ox ^ -SText 

^ ■ a ■ 5^’ 


^9 — ^- Po 


C. = b^. (11) 


An important implication of the scaling relationships for stress, length and dislocation 
density is that any patterns that may emerge from a simulation based on the scaled 
equations automatically fulfill the similitude principle (cf. [TJJ]) - provided that it can 
be demonstrated that the patterns do not depend on initial and boundary conditions 
(i.e.. on system size and deformation history). 

In the scaled variables the kinetic equations read: 

diPi = - 


V„ - (vQ) + vqt. 

( 12 ) 

Vr.(vpt), 

(13) 

V„-('-rQ"’+A®-V„c'), 

(14) 


where the v’elocity v is given by and the quantities ^ and ^ are obtained by 
replacing p,qt aud 5 in Q with their scaled counterparts. The scaled stress is given by 

f = 1 , case ( 1 ), 

f = l + j g{x-x,y-y']^{x'.y')x'y, ca.se(2), (15) 

where the Fourier representation of the kernel Q is given by 


5 (^ 1 , ky) — —Qq 


(16) 


where k = yjk\ + k^ and kx, ky are the components of the w'ave vector in the xy~ 
plane. The dislocation density vector relates to the scaled strain by q = — Vn7, fi^d the 
evolution equation for 7 is 

dn = ^v. (17) 


It is interesting to note that, in the scaled formulation, the model has only two 
uon-dimeusional parameters. The first. ?^,) = aGby/pQ/rext, measures the fraction of the 
externally applied stress that is needed to ov'ercouie short-range dislocation interactions. 
The second parameter, Qq = Gbkpo/r^ measures the relativ’e magnitude of the long- 
range dislocation interactions compared to the external stress. 
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To analyze stability of the model, we envisage a spatially homogeneous reference state 
characterized by the initial values ph = l^qu = 7 /t = with dislocation 

velocity Vh = — T-t where 71 = TioV^ ^^rid with time evolution 


^1% = PhVh , diph = qhVh , diqh = 0 . 


(18) 


We now study the evolution of small fluctuations around this reference state. We write 
these fluctuations as 


7 = + ^7 1 


^ = ph + Spt, 


9t = 9to - 2^n(^7) + 


( 1 ») 


Here, the operator An = Vn * Vn is the 2D Laplacian in the slip plane, and the 
corresponding term in the expression for the curvature density can be envisaged aa 
the curvature of the geometrically necessary dislocations. In linear approximation, we 
And that 

Q“’ = , a'"’ = i(l + Sp,)l. (20) 

The linearized evolution equations for Spi.Sq^, and ^7 are then given by 


di 




Spi 


Sqt 



—Qq 


*5 


1-*/ k J 


0 


2 


0 

i-n 

0 



( 21 } 


Eigenvalues of this matrix are A = A{qt^,7i.QQ. k). We fii'st consider the hypothetical 
case where long-range stresses are appreciable but short-range stresses are absent 
{7{ = 0. ffo > 0). F7om the Routh-H\irw'itz critenon it follows in tliLs case that for 
any values of k. all Eigeuvalues have negative real pail. The system is thus absolutely 
stable and no patterning can occur. More interesting is the opposite case. = 0i ^ >0, 
In this case the Routh-Hui'witz criterion yields the condition 1 > Ti > Tii, where 
= (9 — \/57)/2 as 0,725, for eigenvalues with positive real part to occur, Within 
this unstable regime, the critical eigeuvalue branch is real for all k > 0 , indicating that 
fluctuations on all scales become imdamped while only the homogeneous state = 0 
remains marginally stable. Finally, if both > 0 and 5o > 0, the instability condition 
remains unchanged, however, in this case long-wavelength fluctuations are damped and 
only fluctuations above a critical value iv can become unstable. The regime of unstable 
k vectors is given by 


k > k^iO) 


Gq[ 1Holqt^ sin^+ 

l-i'. 


( 22 ) 


where Hj /2 = (9 ± \/57)/2. Examples showing the critical eigenvalue branch as a 
function of the wave-vector are shown in Fig. 

Since the homogeneous reference state is time dependent with increasing dislocation 
density, one may ask whether instability in the above sense is a sufficient criterion for the 
emergence of heterogeneous dislocation patterns - density fluctiiations which grow more 
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Figure 1; Critical eigenvalue branch as a hmction of k = {kx,ky)\ initial curvature 
= 0.1; left: only short-range interactions {7{ = 0.85.5 = 0), center: long- and short- 
range interactions [Ti = 0.85.5 = 1), right: only long-range interactions (^ = 0,5 = 1}, 


slowly than the statistically homogeneous reference density might in practice hardly 
be visible. While the ultimate verdict about patterning must consider the non-linear 
dynamics of the system and hence be provided by numerical simulation, one can take 
the growth of the reference density into consideration in a linear stability analysis. To 
this end one simply considers, instead of the absolute dislocation density fluctuation 5pi, 
the relative fluctuation Spi/pt,. It turns out that this modification leaves the instability 
condition 1 > “Ho > 'Hi unchanged, but reduces the apparent grcwth rate of the patterns 
and decreases the unstable range of wave vectors lor 5 > 0 as the critical wave vector in 
Eq, ) |2^ is increased by a factor of 2. 

In all cases where \mstable modes occur, the \mstable mode is characterized by 
antiphase fluctuations of 7 and pi (regions of high dislocation density coirespoud to 
regions of low strain and vice versa). The same is true for the flirctnations of 
and pi, hence, the dislocation curvature (and thus the multiplication of dislocations) 
decreases in regions of high dislocation density and increases in regions of low density. 
This is consistent with the well known fact that, in heterogeneous dislocation patterns, 
multiplication occurs by the expansion of segments in the dislocation-depleted regions 
(channels or cell interiors), 

4. Simulation of pattern formation 

The analysis of the previous section can establish conditions for a homogeneous 
dislocation arrangement to become \instable. However, it cannot give access to the 
moiphology of the emerging patterns which depends on the non-linear dynamics of 
the system. Moreover, the instability is in linear approximation not characterized by 
any dominant wavelength, hence, also the characteristic length of the patterns remains 
undetermined by the linear stability analysis. To access pattern morphology and pattern 
scale (if any), we resort to numerical simulation of the model. 
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Fbr the nimierical solution of the scaled CDD evolution equations (]12[>»(|14[>, ( |17| | we 
use a Galerkin finite element scheme together with an implicit time integration. The 
spatial discretization consists of approximately 50 triangiilar elements in each spatial 
direction. In order to properly represent first and second derivatives ^vithin the finite 
element sdieme w use third order Lagrangian shape functions. For reasons of numerical 
stability a small viscous term was added to each of the evolution equations, e.g. in • 12 
we replace the divergence term Vn • { vq ) with * {vq + fV^), where f is chosen 
sufficiently small such that the physical behavior of the equations is not affected but 
numerical oscillations are suppressed. For each time step we update the velocity field 
(i.e. the Taylor term) and then evolve the dislocation fields. One of the benefits of the 
dimensionless scaling is that we can solve systems with widely different parameters based 
on the same finite element mesh: the small but artificial viscous terms for stabilizing 
the numerical sclieme also have the same influence on the residts for widely different 
systems, which improves our control over the numerics. 


4 . 2 . Initial values 

Initial values for the dislocation fields should be consistent with representing a sj'stem 
of curved and connected lines in a coarse-grained manner. To this end we proceed 
in analogy with the coarse graining of a stochastic point process into a superposition 
of Gaussian ’blobs’; we consider a system of discrete dislocation loops of radius r 
and ’smear out’ each loop perpendicular to its line tangent by using a Gaussian 
standard distribution function exp(— 5 ^/ 2 s^) where the standard deviation 

8 characterizes the width of the distribution and f is the distance perpendicular to the 
loop. After having obtained for one loop in this w’ay we determine the corresponding 
dislocation density vector by multiplicatiou of p[ with the radial imit vector, and the 
curvature density by division of by r. Initial values for the fields are then obtained 
by randomly distributing sncli loops and summing \ip the loop fields p[. and q[. An 
example of a resulting initial structure is shown in Fig. ^ An alternative method to 
define initial conditions is to define random scalar fields 7 , and p^ as superpositions 
of randomly located Gaussian ’blobs’ and determine ^ as the gradient of 7 , which 
automatically fulfils the solenoidality condition diva of the Krouer-Nye tensor. To 
make sure that our results do not depend on the construction of initial conditions, both 
methods have been used in the following simulations. 


4‘S. Parameter values 

We consider the system which is defined by the parameters given in table The scaled 
quantities were obtained from the physical values according to Eqs. ([^, (jll]). 
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(a) Pt 



W qt/fh 



0 ) 



luitiaJ Values 


(b) 





Stationary State 


(b) ^ 




Figure 2: Evolution of dimeusionless ODD field variables, for parameters see Table 
ail fields are represented in a projection on the slip plane, instead of gt^e show the more 
intuitive curvature q^ /p^ whicii is the inverse of the local dislocation cur\'atnre radius, Pu 
denotes the edge component of the GND density vector. Surfaces are plotted as elevated 
to tisnallze the fluctuations and are not up to scale. Top: initial field \'alues as obtained 
from a random distribution of 50 dislocation loops, l)ottom: final stationary state; note 
that the scales have changed between top and bottom graphs 
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physical value 

scaled value 

Model parameter 



sv’stcm si2e 

^ =2.00'10-® m 

1 = 17,72 

svstcm height 

/i=0,20'10-® m 

h = 1.77 

Tkylor £attor 

a =0.30 

- 

applied Htress 

=0.17 10® Pa 


total time 

=0.05'10-® s 

fft = 18.97 

Material parameter 



Burgers vector modulus 

6=0,250 10-®m 

6 = 2.27'10-2 

Drag coefficient 

il =5.00'10“^ Pfi 'S 

fl =2.27 ' i0"2 

Shear modulus 

G=128>10® Pa 

C = UG9 

Initial values 



loop radius 

r =0.20 * 10"® m 

r =1.772 

number of loops 

^■ = 50 

- 

average deuslty 

A> =7.85* 10^^ 

Po =1 

btaudaid deviatiou 

.^=25'10-® m 

a =0.222 


Tal)le 2: System and material parameters used for the simulation shown in Fig. 
Dimensionless (scaled) parameters and \'alii€s are indicated by a tilde. 



Figure 3: Evolution of average dislocation quantities in physical units: (a) average 
dislocation density (pi) and number of dislocation loops (b) average velocity (u), 
maximum and niinimimi velocity, and average dislocation flux (p,i') (mean strain rate). 


4.4- Simulation results 

We first look at statistically averaged quantities [volume averages) l)efore we then 
analyze the system in terms of microstructural patterning. 

The mean dislocation density (pt) and loop number are shown in Fig. 15 os 
fmictions of time. We observe that initially the dislocation density quickly increases 
mitil it sat mates at a value that is al)o\it fom times the initial value. At the same 
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Figure 4: Evolution of fluctuations of strain 7 , statistically stored dislocation density 
P.HSD = Pt ~ dislocation velocity i' as ciiaracterized by the respective coefficients 

of variation (COV); the dashed line indicates the point where the instal)ihty condition 
H = 0.725 is first met. 


time the total mimber of dislocation loops Nd stays exactly constant. This is expected 
since loop creation by Prank-Read sources and loop merging by mutual annihilation 
of segments from different loops are not accomited for in the present simple model. 
Accordingly, is ^ conserved quantity - a fact whicli is immediately obvious from Q. 
The fact that stays perfectly constant in our simulation indicates that the numerical 
scheme performs well. 

In order to understand the satiu'ation of the dislocation density, we note that the 
system can become stationary only by reducing the dislocation fiux PtV - and thus the 
strain rate - everywhere to zero. This is possible if either the Taylor stress balances the 
external stress, or if the local dislocation density vanishes. This observation agrees with 
the fact that, in the dislocation-rich regions, the local dislocation density approaclies the 
saturation density given by = [rext/(iiC 6 )]^ = 2.99* 10*^ m“^ where the Taylor stress 
equals the applied external stress. The dynamics of the system is further illustrated in 
the right graph of Fig. which shows the average dislocation flux (the total strain rate) 
and the mean, maximum and minimum dislocation velocity. Initially, we see an increase 
in dislocation flux “which is driven by the increase in dislocation density, This is followed 
by a decrease in dislocation flux caused by the decrease in dislocation velocity due to 
strain hardening. Even in the final state, there remain regions of non-zero dislocation 
velocity which, however, no longer produce strain since they are depleted of dislocations. 

Fig. bottom, shows the CDD field variables in the stationary state; the 
supplementary movie additionally shows the time evolution of the total density. Fig. [J 
shows the concomitant evolution of the fluctiiations of strain, dislocation density and 
dislocation velocity as cliaracterized by the respective coefficients of variation (C'OV, 
defined as the ratio between variance cr and mean of the respective variable). In the 
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initial period (cf. mo^ie), the system is stable since the applied stress is significantly 
higher than the Taylor stress ["H < Tic = Accordingly, fiuctnations of strain or 

dislocation density arc initially smoothed ont as seen from Fig.[^ However, the increase 
in dislocation density driv'es the system towards the critical state which is reached, for 
the simiilated parameters, at a density of about 2.2 x 10*^ in“". This point in time 
is marked by the dashed line in Fig. The ensuing instability is characterized by 
an increase in the fluctuations of strain and dislocation density, and an accelerated 
growth of the fluctuations in dislocation velocity. Once the critical state is readied, 
persistent dislocation density ’blobs’ emerge from the fluctuating density distribution 
while the channels between these blobs become depleted of dislocations. In these 
ciiannels. dislocations move rapidly (the maximum dislocation velocity increases, see 
Fig, 1^ and curved line segments/loops expand and create additional line length which 
is deposited in the blobs. The emergent pattern is characterized by the following features 
(see Fig. bottom): 

(i) Inside the dislocation-rich patdies dislocations are to a significant pari ’statistically 
stored’, i.e. the geometrically necessary dislocation’ (GND) density 1^1 is small. 

(ii) Dislocations that form the boundaries of these dense patches are aligned to each 
other and hence become GNDs, 

(iii) The channels’ between the dislocation rich regions exhibit the largest plastic strain. 
Regions with high density, on the other hand, exhibit a reduced plastic strahi. 

(iv) The curvature is reduced inside the dislocation-dense patches but is high in 
the dislocation depleted ’channels’. The dislocation dense patches act like hard 
inclusions which constrain dislocation motion. Particularly high curvature is 
ol)served in spots where "corners" of the dense patches force the dislocations to 
assume strongly bent configurations, 

(v) The emerging pattern morphology does not depend on the initial values: a smaller 
initial density is compensated by a longer simulation time (more multipMcation due 
to loop expansion) leading to the same final density and morphology. Simulations 
with different type of initial conditions {random fields vs. smeared-out loops) also 
produce the same final morphology’. 

(vi) The emergent pattern morphology does not depend on the numerical discretization 
(as long as the mesh is small enough to be able to represent the patterns shape): 
the pattern l)orders in general do not coincide with the finite element faces, 

(vii) Patterns formed in different simulation runs with the same shear stress exhibit 
approximately the same characteristic length. 

The last obsen’ation is an auspicious feature of all o\ir simulations - regardless of the 
chosen initial values, system size or material parameter. We analyze this behavior in 
the following section in terms of the ’law of similitude’, 
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Patteni evolution under changing stress and the piinciple of similitude 

In Fig. I^a) we show density patterns in simulations where we increase the external 
stress in a step-wise maiiner. such that each stress step is followed by a relaxation step 
until a quasi-stationary configniation is reached before the next stress increment. We 
record the saturation density after each relaxation step. The pattern morphology 
after relaxation remains essentially unchanged diiring the step sequence, however, the 
pattern size decreases with increasing stress while the dislocation density increases. To 
quantify this, we have determined the characteristic length of the patterns. This was 
done by use of the open source image analysis software ImageJ [23|, which allows to 
determine the total number n, area fraction A^ and perimeter Si of each density feature 
i of the density distribution in a semi-automatic way. The characteristic length is then 
obtained as the average value A = (4-4j/S,) of the characteristic sizes of all relevant 
density blobs’. These data are plotted in Fig.j^b) together with the saturation density 
p^'. Fig. I^a) show's examples of the pattern morphology for the largest and smallest 
stresses. Fitting the saturation density data by {fP*'") = [ 7 ^* 1 /({iG6)]^ with a = 0,3 and 



(a) Plot (A) corrcspouds to a shear stress 
rjjjft = 170 MPa, the plot (B) is at r^xt = 
560 MPa. 




(b) Tbe two lines tv'ere fitted with constants a = 0.3 
and D = 9,5. 


Figure 5: Relations between stress Tg^ii patterning length A, and saturation dislocation 
density (pi) (right) and evolution of dislocation density patterns tvith increasing stress 
r,.xt (left). 


the measured pattern wave length by A = DaGb/r^i w'ith D = 9.5 we obtain good 
agreement for ail data; thus the patterns are consistent with the simihtude principle. 
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5. Discussion and Conclusions 

We have formulated a model that desciibes the emergence of heterogeneous dislocation 
patterns that are consistent with the siniilitiide principle. Lineai stability analysis 
has demonstrated that the only requirement for the emergence of such patterns is the 
presence of short-range dislocation interactions as described by the Taylor relationship. 
Long-range internal stresses, on the other hand, are inelevant for the occurrence of 
patterning - although they may influence the pattern morphology by suppressing long 
wavelengths. 

The critical condition formiilnted in Sect, 0.725 < "H < 1. implies a lower and 
an upper boundary to the range of stresses where patterning can occur. If is too 
small (Text < ^r)i neither plastic flow nor dislocation patterning take place. This may 
appear trivial, but it formally demonstrates that the present model is inconsistent with 
the idea of quasi-equihbrium patterning due to energy minimization in the absence 
of externally dri\'en plastic flow. If is too large, i.e, if the deformation rate is 
too high or, more generally speaking, the rate-dependent contribution to the flov,' 
stress becomes too large, then patterning is also suppressed. This observation simply 
means that the flow stress must be predominantly (to at least 72,5 %) attributed to 
dislocation interactions, and not to interactions with the crystal lattice [rate-dependent 
drag component). Qualitatively this is consistent with the fact that heterogeueous 
dislocation patterns are almost generically observed in fee metals, but not in bcc metals 
at low temperatures where the flow stress is controlled by the interaction of dislocations 
with the lattice, 

A mimber of DDD simulations reported in the literature And that, in order to 
produce effective dislocation entanglement, the presence of cross slip is essential (see e.g. 
[T51 fT7|l. The same models show little or no significant increase in dislocation density 
with strain (work hardening) in the absence of cross slip. If. as we perceive it. dislocation 
patterning is essentially a corollary' to work hardening, and DDD simulations require 
cross slip to produce hardening, then cross slip is in these simulations a prerequisite for 
dislocation patterning. We leave it open whether the same is true in real specimens since 
we think the role (or not) of cross slip in hardening, notal)ly in multi-slip conditions, 
needs further investigation. 

It is interesting to compare the critical condition found in the present case with 
the results of Zaiser ^ Sandfeld [TJJ] who carried out a similar analysis for a sv'stem 
of straight parallel dislocations. In that case, the instabihty criterion has a similar 
structure. In fact, the mechanism underlying the instability is in both cases the same: 
The dislocation density evolution equations have the structure of conservation equations. 
Therefore, in regions of reduced dislocation velocity, the dislocation density increases 
and vice versa. Due to strain hardening, a local increase in dislocation density in t\irn 
reduces the local dislocation velocity, leading to instability. We thus find that patterning 
is a direct consequence of the short-range dislocation interactions which control strain 
hardening and of the kinematics of dislocations. No other ingredients are needed. 
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What about dislocatiou multiplication? The present model uatnraJly accounts for 
the increase of dislocation line length due to loop expansion. If we compare the resulting 
instability condition, 0.725 < K < 1, with the counterpart deriv'ed for straight parallel 
dislocations by ZaLser ^ Saudfeld |19| . which can be written as2/3<‘H<l, we find 
that in the latter case the unstable baud is more wide. The difference is due to the fact 
that in the present model which considers curved dislocations, regions of low dislocation 
density develop increased curvature and thus dislocation multiplication localizes in the 
dislocation-depleted regions - a fact which is well known from iii-situ obsen’ations of 
dislocation motion in heterogeneous dislocation patterns. At the same time, dislocation 
multiplication is reduced in dislocation dense regions and this effect tends to counter-act 
the instability. 

On a ’macroscopic’ scale (i.e, the system scale) the average flow stress at saturation 
is slightly lower than in the case of a homogeneous dislocation distribution [Fig, 

This is caused by the much lower flow stress in the nearly depleted ’channels' between 
the blobs and has been noted already by H. Mughral)i in |14) , 

Simulations of the model with increasing stress demonstrate that the patterns 
evolve in a manner that is consistent with the principle of similitude, i.e,. the pattern 
wavelength is proportional to the mean dislocation spacing in the fully developed 
pattern, with a proportionality constant D ^ 10, and decreases in inv'erse proportion 
with the stress. All these featiu'es are consistent with obser\'ations. 

At the same time, it is important to emphasize the very substantial limitations 
of the presented model. With this study we do not aim at realistically reproducing 
dislocation patterns observed in any real deformation experiment. Our aim is instead 
very much on the conceptual level: We want to present a minimal model which can be 
considered the skeleton of any realistic patterning model, while leaving it to future work 
to put in the realistic details. More realistic models will need to account for: 

• Boundniy' conditions typical of real deformation experiments. We chose rather 
artiflcial boundary conditions to represent, for the purpose of analysis, the 
possibility of a situation without long-range internal stresses (case 1). or with only 
long-range stresses (case 2). In reality, of course. l)oth features will always be 
present. 

• In this v.'ork we only considered l)nlk-like behavior. Finite-size effects as e.g. 
ol)served in |2 j3 H] for single crystal pillars will play a role as soon as the pattern 
wavelengths approach the system size. Given that dislocations in onr simulations 
generally have travelled al)out 1 wavelength in distance before they get frozen in 
the patterns it is very likely that the system needs to be significantly bigger than 
the pattern size, and that in small s^'stems dislocation loss through the surfaces 
will inhibit patterning. 

• Though heterogeneous dislocation structures do form even in single slip (i.e. in 
hardening stage I), situations with multiijle active slip s^'stems are much more 
important. Accordingly, one needs to account for latent hardening terms coupling 
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the dislocation densities on different slip systems in the expression for the short- 
range interaction stress, and for a more complex structure of the long-range internal 
stress field. 

• The present model accounts for strain hardening in terms of the dislocation density 
increase that arises from loop expansion. Hcrtt'ever, in present form it does not 
represent a realistic hardening model: the reason is that according to ([^ the 
number of dislocation loops is a conserved quantity. This leads to a parabola¬ 
shaped hardening curve, T( <x y/^, which is not foimd in any real material. The 
reason is simply that in real materials the loop number increases by the action of 
Prank-Read sources, whereas at large strains and/or dislocation densities also the 
merging of loops by annihilation of encountering segments needs to be taken into 
account. All these features, while needed for obtaining realistic hardening curves, 
are however incidental to the patterning mechanism. 

• In a realistic model, additioual effects such as the influence of dislocation 
curvature/line tension on the dislocation velocity may need to be taken into account 
hi an explicit manner. 

As a consequence of these limitations, the blob-like pattern morphology oliserved in the 
present model differs from the morphology of dislocation patterns commonly observed in 
real samples. Other models proposed in the literature (see Sect. produce much nicer 
patterns - but they cannot be shown to be in agreement with the similitude principle 
or to be based upon a realistic description of dislocation interactions, We are therefore 
convinced that the present approach, once the necessary model adjustments have been 
made, can seiwe as a robust framework for future investigation of dislocation patterning 
phenomena in a wide range of deformation settings, 
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